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Abstract 

We examine the solutions of the equations of motion for an expanding Universe, 
taking into account the radiation of the inflaton field energy. We then analyze 
the question of the generality of inflationary solutions in this more general setting 
of a dissipative system. We find a surprisingly rich behavior for the solutions of 
the dynamical system of equations in the presence of dissipational effects. We also 
determine that a value of dissipation as small as ~ 10 -7 H can lead to a smooth 
exit from inflation to radiation. 
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I. Introduction 



It has been by now around 16 years since the idea of inflation [0 was introduced in order 
to solve some of the underlying problems of the standard big bang theory. Since then, 
inflation has become one of the most important paradigms of early Universe cosmology. 
The main idea of inflation is the existence of a period of accelerated expansion of the 
Universe, when its energy density becomes dominated by a potential energy density V(4>) 
of a scalar field <p (the inflaton). The dynamics of inflation depend on the specifics of the 
many models available describing inflation, but the basic mechanisms are based on the 
equation of motion of a (homogeneous) inflaton field, represented by its classical equation 
of motion 

+ 3#0 + V'(0) = O, (1.1) 

(overdot represent derivative with respect to time and V = ^r-) together with those of 
gravity. 

The basic assumption of most inflationary models is the one related with the so called 
slow-roll approximation: at some initial time to, the scalar field has a value far from the 
minimum of the potential V(4>). At this time the scalar field energy density dominates and 
the Universe enters the inflationary era. Sufficient inflation, required to solve the main 
problems of the standard big bang theory, demands that rolls down to the minimum 
sufficiently slowly. It is a common belief that during this stage the particle like matter 
components are all rapidly red-shifted away (as well the temperature, if we initially started 
in a thermal bath) and that the evolution of <fi can be well represented by ( |1 . 1| ) . Later, 
when <fi hits the bottom of the potential and starts oscillating in a time scale shorter than 
the Hubble time, all the energy stored in the oscillating field <fi is then released in the 
form of light particles that may be coupled to. These decay products then thermalize, 
reheating the Universe. This final stage of inflation is the so called reheating stage [for 
an account of new developments in the theory of reheating, see, for instance, pj and 
references therein]. 

More recently, however, Berera and Fang || have shown that these two distinct stages 
of inflation can be compounded if the inflaton dissipates its energy into a thermal bath. 
They have shown that this new scenario is still consistent, providing sufficient e- folding of 
inflation and generation of density perturbations of the right magnitude, and that these 
density perturbations are consequences of thermal fluctuations, as opposite to quantum 
fluctuations, generated during the slow-roll stage, in the supercooled inflationary scenario. 
This new picture for inflation has been denominated warm inflation [|], [5|. 

In this modified description of inflation, we must alter ( |1 . 1| ) to take into account the 
continuous dissipation of the inflaton field in a thermal bath. A phenomenological way 
of doing this is, for example, by introducing a friction-like term f](4>)4> in ( |1 . 1| ) . r?(0) is 
the dissipation coefficient that comes from the damping or "decay" of the scalar field 
interacting with an environment or bath degrees of freedom, constituted by, for example, 
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light fields or field modes that may be coupled to. The form of rj(<j)) depends on 
the details of the interaction terms. The inclusion of a friction term in the equation of 
motion of the inflaton field, to study its evolution in the early stages of inflation, is not 
new. It has long been recognized that dissipative processes, associated with the inflaton 
decay during its evolution, are able to slow down the rolling of the field <fi and, therefore, 
that these processes would be able to support inflation in inflationary models such as 
new or chaotic inflation ||. Recent works on microscopic approaches to nonequilibrium 
dynamics of quantum fields |5], [7], H support the introduction of a friction term of the 
above form in the field equation of motion, especially under certain conditions, such as 
near-equilibrium and when a loop expansion is valid. These studies (see also, ||, [L0|]) have 
also demonstrated the generality of dissipation (as well as fluctuations) in the dynamics 
of fields in interaction with a bath (thermal or quantum). Thus, from the moment that 
we consider that the inflaton is coupled to light fields, which is a necessary condition for 
describing reheating, we must also consider the possible effects of backreactions of these 
fields, which include dissipation through the decay of the inflaton also during the stage of 
inflation. 

In principle we may consider, among other things, that quantities such as sufficient in- 
flation, the magnitude of initial density perturbations, the strength of coupling constants 
of the inflaton field with other fields during the inflationary era, not only constrain, but 
also dictate how important dissipation is during inflation. For instance it is known ||T| 



that small coupling constants are necessary in order that the slow-roll method to be a 
valid approximation, as well as for obtaining the right magnitude of initial density per- 
turbations. Equivalently, density perturbations constrain quantum fluctuations to remain 
small during the slow-roll stage of inflation. On the other hand, in the warm inflation 
picture of ||, |, , thermal fluctuations, related to quite strong inflaton dissipation during 
slow-roll, have been shown to be compatible with generation of initial density perturba- 
tions, with minimal requirements on the inflationary potential, besides that of the slow-roll 
approximation. 

Recent developments in nonequilibrium dynamical processes in the early Universe, 
in particular during reheating, and the warm inflation picture, lead us to consider some 
basic questions concerning dissipation during inflation. Not only a proper microphysical 
approach is necessary to understand the appearance of a non-negligible inflaton dissipation 
at the early stages and/or during inflation, but we also need to investigate how the 
classical inflaton dynamic gets changed in a dissipative environment. In this paper we 
will be concerned with the latter question, since it can shed some light about how specific 
prescriptions for dissipation may change the dynamics in a pre-inflationary era, and then, 
possibly, help to constrain field theoretic model candidates to describe inflation in this 
new setting. A detailed microscopic model will be presented in a forthcoming paper. The 
qualitative analysis of the resulting set of field equations is the first step in understanding 
radiation energy properties of the Universe and in discriminating how general inflation is 
in dissipative systems. Such analysis includes the asymptotic behavior of the solutions 
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near the critical points representing the initial singularity and the basic features of the 
trajectories in the phase space that will be important for the choice of initial conditions 
used in the numerical experiments. We perform this study in the framework of the theory 



of dynamical systems [12|. We then generalize earlier works on this subject [see the Refs. 
[|I3| , which were previously done in terms of the classical set of equations for gravity 
plus that of the inflaton field 0, without dissipation. This study is important both to 
understand the dynamics of inflation in a scenario like warm inflation, as it is also an 
example of a dissipative dynamical system, which may have other important applications 
in early Universe cosmology. 

The paper is organized as follows: In Sec. II we give the basic equations, define 
the region in phase space of interest to our study and discuss some of their properties for 
different kinds of dissipation. In Sec. Ill we present an analysis of the dynamical system of 
equations in the presence of dissipation and determine the relevant asymptotic solutions 
of the system of equations. In Sec. IV we study how dissipation changes the number 
of physically relevant (advantageous) trajectories, i.e., those having, e.g., a long enough 
inflationary stage. We also give limiting values of when dissipation is of relevance in the 
inflationary scenario. Sec. V is devoted to discussions of our results and conclusions. 



II. The Dynamical System 

We start by writing the basic equations that characterize our dynamical system: 

o 8ttG k , 

H 2 = —-( P(t> + p Tad )--i, (2.1) 

2H + 3H 2 + ^- = -8ttG(p + p rad ) , (2.2) 
cr 

<j> + 3H<}> + V'((f>) + 7j(0)0 = , (2.3) 

where H = a J a is the Hubble parameter and G = l/M^, with M p i the Planck mass. 
k — 0, +1, —1 for a flat, closed or open Universe, respectively. pM Ta d) and p^( ra d) are the 
energy density and pressure for (radiation), respectively. We also have the standard 
relations: = |0 2 + V((f>), p^ = \4> 2 — V(4>) and p ra( j = |p ra d- Equations fl2.1|) and 
(|2.2|) are the standard gravity equations derived from the Einstein equations. Eq. ( EOl ) 
is the equation of motion for the (homogeneous) scalar inflaton field 0, with V(4>) the 
effective potential for 0. We have introduced dissipation through a friction term rj(4>)(f) in 
the equation of motion. 

Equation ([O]) deserves some important considerations. First, in a realistic model 



the field dependence of the friction term may not be exactly like the one given in (2.3), 



particularly in highly non-equilibrium situations and in a time dependent background, 



as some recent studies have indicated pi. It has been demonstrated [p], H, however, 
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that 7/(0)0 describes quite well the dissipation of the field in a bath environment in 
near-equilibrium and when a perturbative expansion is valid. In spite of this and without 
a better understanding from a field theoretic point of view of nonequilibrium systems , 
we take for simplicity the friction term in the form given in Q2.31 ). Since our analysis 
is mostly qualitative, we expect that we will not miss much of the important behavior 
of the system given a different prescription for dissipation. Secondly, we take 77 (0) as 
a polynomial in 0. This is the case for standard polynomial interactions among fields, 
at least when a loop expansion in the fields is valid. The field dependence of 7/(0) is 
highly dependent of the type of interaction we consider of the scalar field with the 
fields (or modes) that describe the bath environment. For example, rj(<p) may in some 
approximations be well represented by a constant in the case of linear couplings (like 
4>ip) 0, or be quadratic in for quadratic couplings (like (p 2 ip 2 ) 0. r](<f)) may also 
have a nontrivial dependence on coupling constants and temperature, as is the case at 
high temperatures, which is the situation of interest in warm inflation scenarios. Further 
considerations on 77 (0) will be considered below. Finally, Eq. ( |2.3|) must be viewed as the 
equation obtained by averaging over a noise field, which is directly related to dissipation 
through some generalized fluctuation-dissipation relation. This means we are considering 
the inflaton field as a classical background || and we are looking at the averaged behavior 
of the inflaton field in time. This is a good approximation as far as we look to the inflaton 
dynamics at a sufficiently long period of time when compared to the typical "wavelength" 
of the fluctuations, given by ~ [V]" 1 / 2 . Our scale of time will be basically the duration 
of inflation, which is at least 20 times larger than that value (for e-folds of inflation larger 
than 70), as we will see in Sec. IV. Also, the amplitude of fluctuations, given by the noise 
field two-point correlation function, should not be too large during the inflaton evolution, 
implying, roughly, that our approximation of taking the ensemble averaged equation of 
motion is better the smaller is the dissipation. Keeping in mind these considerations, we 
proceed with the definition of the equations representing our dynamical system. 
We have that p^ and p ra d evolve in time as (p: 

+ 3#0 2 + r7(0)0 2 = , (2.4) 
p rad + AH Pvad - t/(0)0 2 = . (2.5) 



Assuming a flat Universe (k — 0), from (2.1), we can consider 



Prad = 2 -P, = - h 2 - V{<j>) (2.6) 

87TG 07TCt L 



as the first integral of (2.5). Using p^ and p ra d in (|2.2|) , we can write for H the equation: 



H = -2H 2 - (0 2 - 4U(0)) . (2.7) 
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The equations (pl3|), ( |2.7| ), together with = define a three-dimensional (dissipative) 
dynamical system in the phase space of (0, </>, i7). 

For definitiveness, we take the simplest potential V(4>), which describes a massive 
scalar field (f>: V(4>) = \m 2 ^) 2 . The generalization of our results for more complicated 
potentials, including self-interactions terms and symmetry breaking, should not present 
too much difficulty. Next, we take r?(0), as discussed later, in the form = rj n (j) n . 
We can further constrain the form of rj(4>) by noting that the friction coefficient must be 
positive definite due to entropy requirements and the parity invariance of the potential. 
Thus, t] n > and n = or an even integer. We analyze two main cases of possible interest: 
that of n = 0, in which case r] reduces to a constant; and that of n = 2, when rj oc <fi 2 . 
Substituting V and n in ( f2~3l) and (p.7|) and redefining variables t — > t/m, <fi — > ^% nG x , 

6 — > 7=^ y, H — > mz, r/Q — > m/3 and m — > 87rGm/3 2 , we obtain the dynamical system: 

Vo7rG- 



x = y 

y = _ 3yz - x -p nX n y (2.8) 

z = -2z 2 - -y 2 + -x 2 . 

6 y 3 

The physical region in the phase space (x, y, z) is defined by the condition p ra d > 
which, in the units of ( |2.8| ), corresponds to: 

3z " ~ \ y2 ~ ¥ 2 - ' (2 - 9) 

or, expressing the above relation for z, z 2 > |(x 2 + y 2 ), which defines a conical volume in 
(x,y, z), with the physical region p rad > inside the cone. Expanding universes (H > 0) 
correspond to those trajectories for which z > (inside the upper half of the conical 
volume). For our simple quadratic potential there is only one critical point, namely, the 
origin that represents the Minkowski spacetime. In general it is an attractor and the form 
of approach depends upon the dissipation term. In fact we can distinguish two main types 
of approach depending on the kind of dissipation we may have. If n = 0, we can readily 
show from ( |2.8j ) that the trajectories on the physical region approach the origin in the 
direction given by the z-axis. The z-axis itself represents the usual Friedmann model with 
pure radiation and is, in this case, an attractor for the late stages of field evolution. If 
Po < 2 the approach to the z-axis is oscillatory, while for (3q > 2 it is exponential. If n = 2 
the approach is always oscillatory tending to the origin on the surface p ra & = 0. In Figs. 
1(a) to 1(c) we plot the field trajectories by a direct numerical integration of the system 
(|2.8|) for some representative values of initial conditions and dissipative coefficients. 

In the absence of dissipation, the conical surface p ra d = is an invariant manifold. As 
a consequence, all initial conditions chosen on this surface will produce trajectories that 
evolve on it and end at the origin. In particular, the dynamical system for which p ra( i = 
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has been analyzed by Belinskii et. al. [I3[ |H |. The effect of dissipation is to destroy the 
invariant manifold since 



PradU ad=0 = . (2.10) 

In this case the trajectories may no longer lie completely in the surface p ra d = 0, even 
if the initial conditions are taken on it. An important question we can ask is: in which 
way do the curves cross the surface p ra d = ? A simple way to answer this question is, 
for example, by taking the scalar product at p ra d = between the vector normal to the 
outside surface of the cone, n oc — Vp ra d = (x, y, —6z), with the velocity field vector of 
the trajectories on phase space, v = (x, y, z): 

^Uad=o oc -(3 n x n y 2 . (2.11) 

Since we are considering (3 n > and n even, we obtain that n-v\ Ptad =o < 0, which means 
that those trajectories that initially lie inside the cone (p ra d > region) remain inside 
the cone. The only trajectories that can cross the cone lie initially outside it. This is 
a consequence of the requirement of rj((j)) be always positive definite. In particular we 
note that for an odd integer n, we may have ( p. 11] ) positive for x < and we could have 
trajectories leaving the region p rad > to p ra d < 0, in violation of the second law of 
thermodynamics. 

The inflationary region in phase space is defined by 



a 

which, in terms of x, y and z, gives 



- = H + H 2 >0, (2.12) 



z 2 -\y 2 + \x 2 >U. (2.13) 
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Trajectories satisfying (|2.13|) are the ones that are in the inflationary regime. The inter- 



section between the regions defined by ( |2.9|) and (|2.13| ) defines the inflationary physical 
region of interest to us. 



III. Analysis at Infinity and Asymptotic Solutions 

The analysis at infinity provide us with important information concerning the behavior 
of solutions in the region where x 2 + y 2 + z 2 — > oo. Such a region contains information 
about the early stages of the universe just after the quantum domain, and therefore can 
be described by the classical equations. An important issue to be discussed is whether 
inflation plays the role of an attractor for most of solutions no matter the presence of 
dissipation. For our work, the following coordinate transformation is found to be useful: 
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where infinity is represented by the invariant manifold w — 0. The projections of the 
physical region p ra d > as well as the inflationary domain - > on this plane, yield, 
respectively, the regions described as u 2 + v 2 < 6 and 2 u 2 — v 2 > 0. The idea underlying 
the choice of coordinates (|3.1| ) is to compactify infinity, so that a dynamical system in 
these variables describes properly the basic features of the asymptotic behavior of the 
solutions. In Fig. 2, we present a better view of the compactified region given by the cap 
of the conical region. 

The dynamical system (|2.8|) in the coordinates u, v, w is 



(I U fl ( V ! ! 1 ii -L I 

— = uw 2 \2-\ \+vw 2 ^ 

dr \ 6 3 

d D n ( 1) Xl \ " _i I I n 

— = VW 2 -H uw2 +L - (3 n u n vw L 2 (3.2) 

dr \ 6 3 

dw n,, / v 2 

= r+i 2 + -- 
dr \ 6 

where we have changed the time parameter in such a way that dt = w^ +1 dr, which is 
valid for n = 0, 2. The critical points at infinity lie on the plane w = and their nature 
depend on the type of dissipation we are considering. They represent the initial singularity 
since, from Eqs. (3.5)-(3.8) below, we see that the scale factor goes to zero at these points. 
In addition, at the critical points the scalar of curvature R = 6 (H + 2 H 2 ) = ^ (2 u 2 — v 2 ) 
diverges. As stated before two cases will be analyzed, namely, n — 0, for constant 
dissipation, and n = 2, for field dependent dissipation. 



III. a Case n=0 

We adopt the following steps here and for the case n = 2 to analyze the region at infinity 
x 2 + y 2 + z 2 — > oo. First, we determine all critical points on the surface w = that are 
located in the region defined by < 6 < 9* (cf. Fig. 2) (tan6** = a/6), or in terms of the 
new coordinates ( |3.1| ), inside the region u 2 + v 2 < 6. Secondly, in performing the stability 
analysis, we may obtain approximate solutions near those critical points and lying on the 
physical region. Thirdly, the knowledge of unphysical trajectories that lie entirely on the 
surface w = is important since, by continuity arguments, such trajectories determine 
the behavior of the integral curves of the system that lie within the cone but near its 
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After substituting n = into the system fl3.2p and taking w = 0, a two-dimensional 
system arises 
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du ( v 2 

— = u 2+ 

dr \ 6 3 

/ , V 2 

Tt = v h 1 + ir 



(3.3) 



Five critical points are found: two repelling nodes P 2 (u = 0,v = P 2 (u = 0,v = 
— \/6), one saddle P (u = 0, v = 0), and Pi (u = a/6, t> = 0), P[ (u = —y/E, v = 0) which 
are degenerate. In Fig. 3, these points are shown on the plane w = together with the 
unphysical trajectories on it. Note that the projection of the cone p ra d = onto the plane 
w = is an integral curve of the system ( |3.3|) . The regions near P 1( P[ are the projections 
of the inflationary domain 2 u 2 — v 2 > on w = 0. According to the phase portrait shown 
in Fig. 3, the inflationary domain is an attractor for most trajectories. Actually this is 
a general feature of the system, as can be seen from Fig. 4, where we have plotted - 
as a function of time. Initial values were purposefully taken away from the inflationary 
region. Note how the system is driven to the inflationary domain in the different cases of 
dissipation. 

Curves emanate from P 2 and P' 2 to the interior of the cone. The asymptotic behavior 
near these points is 

1 j. , / 3 „ , V M PI ^_,x2_ A) 



H=- ^± lr M pl h(mt) + -^M^^-i, (3.4) 
ot V47T 2v8tt m M p \ 

where Vq is an arbitrary constant and the signs ± indicate which critical point P%, P' 2 
we are considering. The emergence from these points corresponds to the growth of time 
t from the instant t = of the initial singularity. Near these singularities, <fi 2 ^> <p 2 , 
meaning that the kinetic energy of the scalar field is dominant over the potential V(4>) 
producing an effective equation of state = p<j>. Therefore all curves emanating from 
P2, P' 2 are not in the inflationary regime. Note that the influence of dissipation starts to 
be considerable as far as the trajectories evolve away from the critical points. It is also 
possible to integrate the dynamical system near the critical point P , whose asymptotic 
behavior near it is characterized by 

H = ± <f> = ±-^= (3.5) 

where the initial singularity is denoted by t — 0. 

The points Pi, P[ lie on the inflationary region and are attractors of trajectories on 
the invariant manifold w — 0, suggesting that they are also attractors for trajectories on 
the physical region but close to the invariant manifold. It can be shown that there is 
only one trajectory emerging from P x and P[ into the physical region, whose asymptotic 
solution is 
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H=--mt, (f> = ± y=r t, (3.6) 

O Vo 71 

where t increases from — oo, the initial singularity. As expected, this solution is infla- 
tionary, since m 2 2 3> <fi 2 , which produces an effective equation of state p<f, = —p<j>. All 
remaining trajectories that visit a small neighborhood of these points experience an infla- 
tionary phase as described above. It is worth pointing out that the asymptotic solution 
( |3.6| ) corresponds to the direction of approach given by v = — w, which is equivalent to 
no dissipation (see analysis carried out by Belinskii et al ||14|| ). The influence of dissipation, 
from a purely dynamical system analysis, appears as a second order effect changing the 
former relation to v = — ^ w + ^ /3 w 2 . 

IH.b Case n=2 

The structure of the phase space for n = 2 is quite different from the previous case. After 
substituting n = 2 into the system (|3.2[) , we obtain a two-dimensional dynamical system 
whose phase portrait is shown in Fig. 5. There are two lines of degenerate critical points, 
denoted by 1% (— oo < u < oo,v = w = 0) and l 2 (— oo < v < oo, u = w = 0). Here 
the projection of the cone surface p ra d = 0, represented by the dashed circle on Fig. 5, 
is not an integral curve of the system at infinity. By continuity arguments, taking into 
account the unphysical trajectories on w = 0, we infer that the trajectories within the 
physical region, but near the plane w — 0, originate outside the region p rad = 0, eventually 
outside the quantum domain, i.e., where the classical description is valid. Therefore, they 
must be ruled out of our analysis. Nonetheless, some critical points on the lines li, I2 
may be source of trajectories in the physical region. In fact, using standard methods to 
analyze sets of degenerated critical points ||12]| , we come to the conclusion that the points 
at infinity P (u = v = 0), P x (u = VQ, v = 0), P{ (u = -VQ, v = 0), P 2 (u = 0, v = \/6), 

are the relevant ones and belong to a distinct nature from the 
remaining point, as we are going to describe. From the points Pq, P 2 and P' 2 three- 
dimensional pencils of trajectories emerge, whose asymptotic behavior near the points Pq 
and P 2 , P' 2 , are given, respectively, by 

jr = l, ^ = ±^=, (3.7) 

E = YV <P = ±^ M p^( mt )- (3-8) 

In the above solutions the initial singularity is characterized by t = 0, or equivalently the 
points Pq and P 2 , P 2 , and the emergence from the singularity corresponds to increase of 
time from this value. We remark that the solutions are not inflationary, but eventually 
evolve to this stage. 
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The most interesting class of solutions emanate from the points P\, P[ on the line 
l\ (see Fig. 5). It can be shown that a two-dimensional pencil of paths emerges from 
these points along the same plane of the line l\ and orthogonal to the plane w — 0. The 
asymptotic behavior in this region is characterized by v = 0, which implies that the form 
of the solutions is the same as found in (|3.6|). 



IV. Advantageous Trajectories and Radiation in 
Inflation with Dissipation 

We now have a closer look at those trajectories in phase space that may lead to a long 
enough inflationary stage and let us try to quantify these trajectories in a way similar to 
that in [Tj|. During an inflationary period given by At = tf — ti, the scale factor changes 
by 



exp Qf 1 H{t)dt^j . (4.1) 



a(t f ) < n, 

a(ti) 

Let us determine the above ratio in terms of the initial and final values for the scalar field, 
4>i and 4>f, respectively. 

First note that the inflationary epoch can be characterized by those trajectories that 
satisfy 

m0>|0|. (4.2) 

We consider the slow-roll approximation which, in the presence of dissipation, can be 
written as 

m ^ (4.3) 



3H + 7]((j)) ' 

We also assume that ^> p ra d (see [[J and ||). Then, from ( |2.1| ) and ((4.2[), we can write 



Using ( |4.4| ) and ( |4.3| ) in ( f4.1| ), we get the explicit expression: 

2tt , 2 2 x /47 [*t 7}(<f>) 



a{t f ) 

—f- ~ exp 

a(ti) 



M p 2 V 3 7 0I mM pl 



(4.5) 



For r^(0) = we recover the usual expression for the scale factor ratio |], Taking 
rj((f)) = r] n (j) n , n = 0, 2, we finally obtain 
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—f- ~ exp 
a(ti) 



AfJ v 1 t;j ' V 3 mM p i n + 1 



(4.6) 



The term in the above exponential gives the total number, N, of e-folds during inflation. 
We have numerically checked the validity of the above equation. The physically interesting 
trajectories are as usual those with N>70 [|TJ. 

With no dissipation, r\ n = 0, N ~ 70 can be accomplished with fa « 3.4M pl and 
</>/ ~ 0.2M p i Q. In the presence of dissipation, for a fixed N, dissipation decreases the 
values of both fa and <pt. For the values of fa and <f)f, which N ~ 70 with rj n = 0, we see 



from (4J5) that dissipation begins to change considerably N, say by more than 10%, for 
?jo>m, or r72>0.25m/Mp! (in terms of the dimensionless factors f3 n , /3 >1 and /3 2 >0.01, 
respectively). In both cases we have an initial dissipation of order j](fa)>m. Since H m m, 
we must have H<r]{fa) in order for dissipation to begin changing the inflationary phase 
considerably. 

The authors in [|l3|, [ljj have characterized the measure of disadvantageous trajectories 
(N < 70) as ~ m/M p i, obtained by the ratio of the trajectories beginning at \<p\ < fa, 
in relation to the complete quantum boundary, of length 2TTM^ l /m. Since for fixed N 
the effect of dissipation is to reduce fa, we find that in the presence of dissipation m/M p i 
(~ 10 -7 ) is the upper bound for the measure of disadvantageous trajectories. Dissipation 
turns inflation even more general. 

We can also evaluate the energy of the inflaton field transferred to radiation during 
the inflationary stage, from, for example, Eq. (|2.5|) . In Fig. 6 we plot p ra d against 
time for four different values of f3 n . The mark indicates when the trajectories leave the 
inflationary region (a/a becomes negative). We find that there is a value of (3 n for which 
radiation is maximal, at the exit of the inflationary stage, given by /3q, /3 2 ~ 1 (rjo ~ m 
and r)2 ~ 25m/Mpj). These results are fairly independent of the initial conditions. Also, 
we can find an ideal range for f3 n such that we exit inflation smoothly to the radiation era 
with the Universe sufficiently hot. At thermalization 

Prad = , (4.7) 

where g* (~ 10 2 — 10 3 ) is the effective number of degrees of freedom. If we impose that 
at the end of inflation we should have at least T ~ m (m ~ 5.10~ 7 M p i for the quadratic 
potential model Q), we have numerically obtained, for both j3 and /3 2 , the approximate 
bound: 



l(T 7 </3„<800 . (4.8) 

Too high values of j3 n lead to a longer inflationary period. Most of the inflaton's energy 
is transferred to radiation which is red-shifted away before the end of inflation. There is 
no possibility of a reheating phase for these cases. On the other hand, very small values 
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of (3 n can provide a smooth exit from inflation to radiation and the inflaton field energy 
density can still provide an additional reheating phase by the usual mechanisms^. The 
maximum temperature reached by the thermal bath, at the exit of the inflationary stage 
in our model, is T « 10 15 GeV, for /3 , /3 2 ~ 1. 



V. Discussions and Conclusions 

We have examined a dynamical system describing inflaton dynamics in the presence of 
dissipation. Dissipation is taken as being present throughout the inflaton's evolution and 
we have presented the changes in the phase diagram due to the effect of dissipation. Our 
results allow us to draw a series of important conclusions concerning the dynamical system 
and the physical system itself. 

We have shown that the continuous production of radiation due to dissipation changes 
effectively the inflaton dynamical system. We have distinguished two main consequences 
for the field trajectories in phase space due to the presence of dissipation: (a) In field 
independent dissipation (for (3q > 2) the trajectories tend towards the z-axis at the final 
stages of evolution. In this regime the inflaton's motion is overdamped. The stronger 
the dissipation (/3q 2) the faster the field trajectories are attracted to the z-axis; (b) 
In field dependent dissipation the motion is initially damped but the field trajectories 
always tend to the cone surface p ra( j = 0. The stronger the dissipation the slower the field 
trajectories tend to this surface. 

From the physical system side we can also reach some interesting conclusions which 
may have important consequences for inflationary dynamics. We have seen that dis- 
sipation associated with inflaton decay not only supports inflation but also draws all 
inflationary field trajectories into a longer stage of inflation. This is true for both kinds 
of dissipation we have considered. In terms of a measure of the amount of trajectories 
possessing the required inflationary period, a similar analysis to the one given in |L3|, |14| 
allows us to reach the conclusion that inflation is even more general in the presence of 
dissipation. 

Our results also allow us to reach some conclusions concerning warm inflation scenarios 
or non-isentropic inflation, in particular concerning the problem of reheating. In non- 
isentropic inflation, as opposed to standard inflation, there is the possibility of no reheating 
phase, due to the continuous production of entropy and heat, leading to a smooth exit 
from inflation to the radiation era || We then expect that in these scenarios there 
is a strong damping of the inflaton field's motion in the direction of the minimum of the 
inflaton's potential. In this paper we have quantified how strong the dissipation must be, 
when compared to the friction due to expansion, in order to achieve favorable scenarios of 
radiation at the exit from the inflationary stage. We have shown that no reheating is only 



*At least for /3 < 2, when we are in the underdamped regime (see Sec. II) and for any value of /?2 



smaller than the upper value in (4.8) 
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achieved in the overdamped regime (/3q > 2) for constant dissipation, when there is no 
oscillatory motion of the inflaton around the potential minimum, or for very high values 
of /?2) when all the inflaton's energy is transfered to radiation and this is then red-shifted 
away due to the longer expansion. 

A further question we may ask is about a possible realization of chaotic behavior 
in Cosmology. According with some previous in the realm of isotropic and anisotropic 



Bianchi IX cosmologies||16||, the Universe can experience two possible outcomes: collapse 
into the big crunch after an initial stage of expansion and expansion into the inflation- 
ary regime. In this way, chaotic behavior is associated with an indetermination of the 
final state (collapse or expansion) of the Universe once a set of initial conditions in an 
infinitesimal region is chosen. We may say that the boundaries of collapse and expansion 
are mixed. In our caseQ all orbits initially in the physical region p ra( i > are unavoidably 
attracted to the state of unlimited expansion represented by the origin of the phase space 
(cf. Fig. 1). The orbits can not randomically explore the region inside the cone p ra( i = 0, 
since the volume of the phase space is not conserved. Therefore, we can conclude that 
chaotic behavior does not take place for those orbits lying in the physical region in the 
cases studied previously, i. e., for n — 0, 2. Such conclusions are also valid for the case 



analyzed by Belinskii and Khalatnikov |T^] in the realm of anisotropic models. 

Finally one of the most important results of this paper was to show that even a value 
of dissipation as small as rj(<j)) ~ lCT 7 m (compare with H ~ m) can still lead to very 
important consequences for inflationary dynamics, such as providing a smooth exit from 
inflation to radiation, and it also shows that the interaction of the inflaton field with other 
degrees of freedom, manifested in the form of dissipation, cannot in general be neglected. 
We are currently working on a detailed microscopic model for dissipation in inflation, to 
be presented in a forthcoming paper. 
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Figure Captions 



Figure 1: The field trajectories in phase space (projected in the plane x, z) for the cases: 
(a) no dissipation; (b) constant dissipation (fi = 0.3 for dashed line, /3 = 0.7 for full line) 
and (c) field dependent dissipation (/3 2 = 10 for full line and (3 2 = 90 for dashed line). 
Initial conditions for both cases are (x ,y ,z ) = (v6, 10~ 3 , 2.5). Dotted lines represent 
the lines where p ra d = 0. 

Figure 2: View of the physical region of the phase space. The infinity is compactified 
onto the plane w = 0. 

Figure 3: Phase portrait of the plane w = for the case of constant dissipation {n = 0). 
The inflationary domain is represented as the projection into the plane w = 0. 



Figure 4: Plot of - against time, in the units of ( |2.8|) , for the cases of: no dissipation 
(full line); constant dissipation (/? = 10), dashed line; and field dependent dissipation 
(/3 2 = 10/6), dotted line. Initial conditions for both cases are (x , y , z ) = (V6, 10~ 3 ,2.5). 

Figure 5: Integral curves at infinity for the case of field dependent dissipation. The 
dashed circle represents the projection of the cone p ra d = 0. The projection of the infla- 
tionary domain, into the plane w = 0, is the same as in Fig. 3. 

Figure 6: p ra d for different values of dissipation. The mark indicates the exit of the 
inflationary stage. Initial conditions are (0j,0j,ifj) = (3.4M p i,0,3.4 Jym) (Pradi = 0). 
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FIGURE 4 
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FIGURE 6 
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